Separation of Gravitational- Wave and Cosmic-Shear Contributions to Cosmic 

Microwave Background Polarization 
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Inflationary gravitational waves (GW) contribute to the curl component in the polarization of 
the CMB. Cosmic shear — gravitational lensing of the CMB — converts a fraction of the dominant 
gradient polarization to the curl component. Higher-order correlations can be used to map the 
cosmic shear and subtract this contribution to the curl. Arcminute resolution will be required to 
pursue GW amplitudes smaller than those accessible by Planck. The finite cutoff in CMB power 
at small scales leads to a minimum detectable GW amplitude corresponding to an inflation energy 
near I0 15 GeV. 



Observation of acoustic oscillations in the temperature 
anisotropics of the cosmic microwave background (CMB; 
|EJ) strongly suggests an inflationary origin for primor- 
dial perturbations 0. It has been argued that a new 
smoking-gun signature for inflation would be detection of 
stochastic background of inflationary gravitational waves 
(IGWs) . These IGWs produce a distinct signature in 
the CMB in the form of a contribution to the curl, or 
magnetic-like, component of the polarization 0). Since 
there is no scalar, or density-perturbation, contribution 
to these curl modes, curl polarization was considered to 
be a direct probe of IGWs. 

There is, however, another source of a curl compo- 
nent. Cosmic shear (CS) — weak gravitational lensing 
of the CMB due to large-scale structure along the line 
of sight — results in a fractional conversion of the gradi- 
ent mode from density perturbations to the curl compo- 
nent ||. The amplitude of the IGW background varies 
quadratically with the energy scale E- ln f[ of inflation, and 
so the prospects for detection also depend on this energy 
scale. In the absence of CS, the smallest detectable IGW 
background scales simply with the sensitivity of the CMB 
experiment — as the instrumental sensitivity is improved, 
smaller values of E\ n & become accessible |||| • More real- 
istically, however, the CS-induced curl introduces a noise 
from which IGWs must be distinguished. If the IGW 
amplitude (or £i n a) is sufficiently large, the CS-induced 
curl will be no problem. However, as E- m fi is reduced, 
the IGW signal becomes smaller and will get lost in the 
CS-induced noise. This confusion leads to a minimum 
detectable IGW amplitude fflj. 

In addition to producing a curl component, CS also 
introduces distinct higher-order correlations in the CMB 
temperature pattern. Roughly speaking, lensing can 
stretch the image of the CMB on a small patch of sky 
and thus lead to something akin to anisotropic correla- 
tions on that patch of sky, even though the CMB pattern 
at the surface of last scatter had isotropic correlations. 



By mapping these effects, the CS can be mapped as a 
function of position on the sky Q. The observed CMB 
polarization can then be corrected for these lensing de- 
flections to reconstruct the intrinsic CMB polarization at 
the surface of last scatter (in which the only curl com- 
ponent would be that due to IGWs) . In this Letter we 
evaluate how well this subtraction can be accomplished 
and study the impact of CS on experimental strategies 
for detection of IGWs. 

To begin, we review the determination of the smallest 
detectable IGW amplitude in the absence of CS. Follow- 
ing Ref . H , we consider a CMB-polarization experiment 
of some given instrumental sensitivity quantified by the 
noise-equivalent temperature (NET) s, angular resolu- 
duration i yr in years, and fraction of sky 
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covered / s k y . We then make the null hypothesis of no 
IGWs and determine the largest IGW amplitude T, de- 
fined as T = 9.2V/m,p l , where V = E* nfi is the inflaton- 
potential height, that would be consistent at the la level 
with the null detection. We then obtain the smallest de- 
tectable IGW amplitude <jt from 
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where (j® B,GW is the IGW contribution to the curl power 
spectrum, and 
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is the standard error with which each multipole mo- 
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be determined. Here, w = 
is the variance per unit area 



on the sky for polarization observations when i p j X is the 
time spent on each of iV pix pixels with detectors of NET s, 

1 /2 

and 9 ~ 203 / / is roughly the width (in degrees) of the 
survey. In restricting the sum to I > 180/0, we have as- 
sumed that no information from modes with wavelengths 
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FIG. 1. Minimum inflation potential observable at la as a 
function of survey width for a one-year experiment. The left 
panel shows an experiment with NET s = 25 /iK y^sec. The 
solid curve shows results assuming no CS while the dashed 
curve shows results including the effects of an unsubtracted 
CS; we take #fwhm = 5' in these two cases. The dotted curves 
assume the CS is subtracted with (3fwhm = 10' (upper curve) 
and 5' (lower curve). Since the dotted curves are close to the 
dashed curve, it shows that these higher-order correlations 
will not be significantly useful in reconstructing the primor- 
dial curl for an experiment similar to Planck's sensitivity and 
resolution. The right panel shows results for hypothetical im- 
proved experiments. The dotted curves shows results with 
CS subtracted and assuming s = 1 /iK y^sec, Sfwhm = 5', 
2', and 1' (from top to bottom). The solid curve assumes 
0fwhm = l' and s = 1 fj,K y^sec, and no CS, while the dashed 
curve treats CS as an additional noise. The long-dash curve 
assumes CS subtraction with no instrumental noise (s = 0). 



larger than the survey size can be obtained; in fact, some 
information can be obtained, and our results should thus 
be viewed as conservative |7|. 

The second term in Eq. (^) is due to instrumental 
noise, and the first is due to cosmic variance. In the ab- 
sence of CS, and for the null hypothesis of no IGWs, we 
set Cf B — 0, and the results for the smallest detectable 
IGW amplitude are shown as the solid curves in Fig. |l| for 
an experiment with detectors of comparable sensitivity to 
Planck's (left) and a hypothetical experiment (right) with 
better sensitivity. The smallest detectable IGW ampli- 
tude T scales as s 2 t~ 1 . For large survey widths, it scales 
as 9, but at survey widths smaller than ~ 5° it increases 
because information from the larger-angle modes in the 
IGW-induced curl power spectrum is lost (cf. the IGW 
power spectrum in Fig. 0). 
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FIG. 2. Contributions to the CMB polarization power 
spectra. The long-dashed curve shows the dominant polariza- 
tion signal in the gradient component due to scalar perturba- 
tions. The solid line shows the maximum allowed curl polar- 
ization signal from the gravitational-wave background, which 
will be smaller if the inflationary energy scale is smaller than 
the maximum value allowed by COBE of 3.47 x 10 16 GeV. The 
dashed curve shows the power spectrum of the curl component 
of the polarization due to CS. The dotted curve is the CS con- 
tribution to the curl component that comes from structures 
out to a redshift of 1; this is the level at which low-redshift 
lensing surveys can be used to separate the CS-induced polar- 
ization from the IGW signal. The dot-dashed line is the resid- 
ual when lensing contribution is separated with a no-noise 
experiment and 80% sky coverage. 



It is now easy to see how inclusion of CS affects these 
results. As discussed above, lensing of the gradient po- 
larization at the surface of last scatter due to density 
perturbations leads to a CMB curl component with a 
power spectrum, 
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li here and throughout, Cf E is the power 
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spectrum of the gradient component of polarization and 
Cf^ is the power spectrum of the projected lensing po- 
tential |l(J . The latter is defined in terms of the potential 
fluctuations, <E>, along the line of sight such that 
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where r is the comoving radial distance, or conformal 
look-back time, with ro at the last scattering surface, and 
dA(r) is the comoving angular diameter distance. The 
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CS-induced curl power spectrum is shown as the dashed 
curve in Fig. |[ 

By the time these measurements are made, the cosmo- 
logical parameters that determine this lensed curl power 
spectrum should be sufficiently well determined that this 
power spectrum can be predicted with some confidence. 
In that case, the CS-induced curl component can be 
treated simply as a well-understood noise for the IGW 
background. The smallest detectable IGW amplitude 
can then be calculated as above, but now inserting the 
lensed power spectrum, Eq. in Eq. (||) for Cf B . 

The results are shown as the short-dash curves in Fig. 
|l|. When lensing is included, the results no longer scale 
simply with / s k y , s, or t yr , as there is now a trade-off be- 
tween the instrumental-noise and CS-noise terms in Eq. 
(Q). The left panel shows that the IGW sensitivity for 
an experiment with NET similar to Planck's should not 
be affected by CS. This is because the IGW amplitudes 
that could be detectable by such experiments are still 
relatively large compared with the expected CS signal, 
especially at the larger angles that will be best accessed 
by Planck. However, CS will affect the ability of exper- 
iments more sensitive than Planck to detect unambigu- 
ously IGWs, as shown in the right panel. CS also shifts 
the preferred survey region to larger areas, as the IGW 
power spectrum peaks at larger angles than the CS power 
spectrum (cf. Fig. ||). Finally, note that if the CS curl 
is treated as an unsubtracted noise, it leads, assuming a 
no-noise polarization map, to a smallest detectable IGW 
amplitude, corresponding to an inflaton-potential height, 
V 1 / 4 - 4 x 10 15 GeV. 

Now we arrive at the main point of this Letter, i.e., how 
well can the CS-induced curl be subtracted by mapping 
the CS as a function of position on the sky? One possibil- 
ity is that the primordial polarization pattern might be 
reconstructed from that observed by using CS maps ob- 
tained with correlations of galaxy ellipticities fll] . How- 
ever, the source galaxies for these CS surveys are at red- 
shifts z ~ 1, while only a small fraction of the CS-induccd 
curl comes from these redshifts, as indicated in Fig. ||. 
An alternative possibility is to use higher-order correla- 
tions in the CMB M to map the CS-induced curl all the 
way back to the surface of last scatter. 

CS modifies the temperature and polarization pattern, 
giving rise to anisotropic correlations on small scales 
where the image of the CMB surface of last scatter is 
sheared by weak lensing. According to Ref. [I4I , the 
quantity, V • [T(n)VT(n)], provides the best quadratic 
estimator, given a temperature map, of the deflection 
angle at position n on the sky. In Fourier space, we can 
write this quadratic estimator for the deflection angle as 

Ni f d 2 h , r CMB , 1 1 r,CMK\ T (h)T(h) 
( ) = T J (2^ ( ' 1 h +1 * 12 ^ > 2Cjf Cjf ' 

where C 1 ; is the unlensed CMB power spectrum and 



c tot = ^lcnscd-CMB + gl a «r* includes all contri . 

butions to the CMB temperature power spectrum. The 
ensemble average over CMB realizations, (q:(1)cmb)i is 
equal to the deflection angle, l(j>(l), when 




It can also be shown |l2[] that Ni is the noise power spec- 
trum associated with the reconstructed deflection angle 
power spectrum, 

(a(l)d(l')) = (2^) 2 <5 D (1 + 1') (} 2 Ct* + Ni) . (7) 

Here the ensemble average is taken independently over 
realizations of both the CMB and the intervening large- 
scale structure. In addition to these temperature esti- 
mators for the deflection angle, we also use analogous 
ones constructed from the polarization, as discussed in 
Ref. JTJ, although we do not reproduce those formulas 
here. The total noise in the estimator for the deflection 
angle can then be constructed by summing the inverses 
of the individual noise contributions. We thus determine 
the variance with which each Fourier mode of <p can be 
reconstructed. 

With the deflection angle obtained this way as a func- 
tion of position on the sky, the polarization at the CMB 
surface of last scatter can be reconstructed (details to be 
presented elsewhere In the ideal case, there would 

be no error in the CS reconstruction leading to no resid- 
ual lensing- induced curl component. Realistically, how- 
ever, there will be some error in the CS reconstruction, 
from measurement error and also from cosmic variance, 
with a noise power spectrum, c^' 1101130 , given by Ni/l 2 . 
Since cf "scd-CMB > c cmb along the dampmg tail of 

the CMB power spectrum, Ni converges to a finite value 
such that beyond a certain small angular scale, additional 
information from CMB provides no further information 
on CS. This leads to a finite limit on lensing extraction 
and, subsequently, a limit for the amplitude of IGW con- 
tribution that can be separated from CS. 

The lensing reconstruction from CMB data only al- 
lows the extraction of Cf^ to a multipole of < 1000 [fl3|| , 
but there is substantial contribution to the CS-induced 
curl component from lensing at smaller angular scales. 
We thus replace Ni/l 2 by C, when the former exceeds 
the latter at large I. This provides an estimate to the 
noise expected in the reconstructed curl component that 
follows from implementing a filtering scheme where high- 
frequency noise in the CS reconstruction is removed to 
the level of the expected CS signal. The dot-dash curve 
in Fig. |^ shows the residual CS-induced curl component 
(5i)hat remains after subtraction. 

We can now anticipate the smallest IGW amplitude 
detectable by a CS-corrected polarization map by simply 
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using this residual noise power spectrum in Eq. (||). Re- 
sults are shown as dotted curves in Fig. |. The left panel 
shows results as a function of survey size for an experi- 
ment with NET similar to Planck, while the right-hand 
panel shows results for experiments with better sensi- 
tivity and resolution. Since the dotted curves are just 
below the dashed curve in the left-hand panel of Fig. [I], 
we learn that Planck's sensitivity will not be sufficient to 
warrant an effort to reconstruct the primordial curl and 
we would do just as well to simply treat the CS-induced 
curl as a noise component of known amplitude. We can 
expect to improve the discovery reach for IGWs by in- 
creasing the sensitivity and resolution. The right-hand 
panel of Fig. [I] shows results for a hypothetical experi- 
ment with s = 1 /zK y/sec and angular resolutions of 5', 
2' and V . We now see there is a significant difference be- 
tween the dashed curve and the dotted lines suggesting 
the increasing improvement with increasing resolution. 

Finally, to indicate the ultimate limits of this class of 
experiments, the long-dash curve in Fig. [l] shows the re- 
sults assuming perfect detectors (i.e., s = 0). If there 
were no CS-induced curl, then we would have sensitiv- 
ity to an arbitrarily small IGW amplitude, but the exis- 
tence of a CS-induced curl provides an ultimate limit of 
yi/4 ~ 4 x 10 15 GeV, as discussed above. Correction for 
the effects of CS with the CS map inferred from higher- 
order correlations would allow us to access lower IGW 
amplitudes. If there is no instrumental-noise limitation, 
the sensitivity to an IGW signal is maximized by cov- 
ering as much sky as possible, and the lowest accessible 
inflaton potential, ~ 10 15 GeV, is obtained with a nearly 
all- sky experiment. 

To conclude, we have studied the IGW amplitudes ac- 
cessible by mapping the curl component of the CMB po- 
larization, taking into account the effects of a CS-induced 
curl that is either modeled as an unsubtracted noise or 
subtracted with a CS map obtained with higher-order 
correlations. We find that the CS reconstruction is un- 
likely to improve the IGW discovery reach of Planck. To 
go beyond Planck, however, a CS map will need to be 
constructed with temperature and polarization maps of 
higher sensitivity and resolution than Planck. An ul- 
timate limit of roughly I/ 1 / 4 ~ 10 15 GeV to the de- 
tectable IGW amplitude using the techniques considered 
here comes from the existence of finite CMB power on 
small angular scales. There are several possible ways 
this lower limit may be improved upon. First of all, we 
have used only the lowest-order temperature-polarization 
correlations to reconstruct the CS. The inclusion of the 
complete temperature-polarization four-point correlation 
functions and higher-order correlations may possibly im- 
prove the CS reconstruction. We have neglected reion- 
ization in our analysis. Reionization will boost the large- 
angle CMB polarization [H improving the detectability 
of IGWs. Another improvement in this limit may be 
achieved by relaxing our assumption that power in the 



CMB drops exponentially at small scales. If the excess 
small-scale power recently detected by CBI []l6| comes 
from high redshifts, then there will be more small-scale 
coherence patches with which to reconstruct the CS. In 
this case, it is imaginable that a far more precise CS map 
can be reconstructed, but this might require even better 
angular resolution and sensitivity. 

During the preparation of the paper, we learned of 
other very recently completed work by Knox and Song 
]l7t that performs a very similar calculation and reaches 
similar conclusions. This work was supported in part by 
NSF AST-0096023, NASA NAG5-8506, and DoE DE- 
FG03-92-ER40701. Kesden acknowledges the support of 
an NSF Graduate Fellowship and AC acknowledges sup- 
port from the Sherman Fairchild Foundation. 
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